Dark Energy and the Return of the Phoenix Universe 



Jean-Luc Lehners 1 and Paul J. Steinhardt 1,2 
1 Princeton Center for Theoretical Science, Princeton University, Princeton, NJ 08544 USA 
2 Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544, USA 

In cyclic universe models based on a single scalar field (e.g., the radion determining the distance 
between branes in M-theory), virtually the entire universe makes it through the ekpyrotic smoothing 
and flattening phase, bounces, and enters a new epoch of expansion and cooling. This stable 
evolution cannot occur, however, if scale-invariant curvature perturbations are produced by the 
entropic mechanism because it requires two scalar fields (e.g., the radion and the Calabi-Yau dilaton) 
evolving along an unstable classical trajectory. In fact, we show here that an overwhelming fraction 
of the universe fails to make it through the ekpyrotic phase; nevertheless, a sufficient volume survives 
and cycling continues forever provided the dark energy phase of the cycle lasts long enough, of order 
a trillion years. Two consequences are a new role for dark energy and a global structure of the 
universe radically different from that of eternal inflation. 



Lemaitre [l| invoked the mythical "phoenix" to refer 
to an oscillatory cosmology in which the universe under- 
goes regular periods of expansion and contraction with 
a big crunch/big bang transition in between. His model 
was a closed universe in which the reversal from expan- 
sion to contraction is caused by an overdensity of mat- 
ter and radiation. Although this model is ruled out by 
WMAP measurements Q indicating a flat universe, there 
remains the more recent proposal of a cyclic model of 
the universe [1, 0] (see [f| for a recent review) that also 
undergoes regular cycles of evolution beginning with a 
big bang and ending in a big crunch. Unlike Lcmaitrc's 
model, the cyclic universe is kept smooth and spatially 
flat by a period of ekpyrotic contraction [6J, |7[ (slow con- 
traction with equation of state w 3> 1) that precedes each 
big crunch. The expansion, contraction and smoothing 
are governed by a scalar field 4>\ that evolves along a po- 
tential V(4>\). In a version inspired by M theory Q, the 
cycles correspond to the regular collision and rebound of 
two orbifold planes, and (f>i is the modulus field whose 
value determines the distance between orbifold planes. 
As in Lemaitre's oscillatory model, the notion is that 
virtually the entire universe makes it through the big 
crunch/big bang transition and continues into the next 
cycle. 

In this paper, we show how the generation of scale- 
invariant density perturbations may lead to a variation 
of the cyclic model that may be more fittingly called a 
"phoenix universe," in the sense that most of the uni- 
verse is turned to "ash" at the end of each cycle and 
the universe is reborn from a comparatively tiny surviv- 
ing seed. The phoenix effect is an unintended byproduct 
of generating curvature perturbations via the entropic 
mechanism |9| . The entropic mechanism was introduced 
because an ekpyrotic contraction phase with only a sin- 
gle scalar field produces scale-invariant time-delay per- 
turbations, but not scale-invariant curvature perturba- 
tions Although there are various proposals for cir- 
cumventing this problem , the entropic mechanism 
is currently the best understood and can be fully an- 
alyzed in four dimensional field theory without invok- 
ing extra dimensions 0, [13, [H, [HI- The mechanism 



relies on having two scalar fields (e.g., the radion and the 
Calabi-Yau dilaton [9() each obtain scale- invariant fluc- 
tuations by evolving down a steep potential V(<f>i, 4>2) 
during the ekpyrotic phase, with one linear combination 
forming a spectrum of time-delay fluctuations and the 
second a spectrum of entropic perturbations. The latter 
are converted to curvature fluctuations if the two-field 
trajectory bends or bounces as the universe approaches 
the big crunch, which occurs generically. 

An essential feature of the entropic mechanism is that 
the classical ekpyrotic trajectory a long V(d>u §2) is un- 
stable to transverse fluctuations conse- 
quence is that, without adding some mechanism to force 
the universe to begin very close to the classical track 
[HI], an overwhelming fraction of the universe fails to 
make it all the down the classical trajectory simply due to 
quantum fluctuations. This fraction is transformed into 
highly inhomogeneous remnants and black holes that do 
not cycle or grow in the post-big bang phase. However, as 
shown below, something curious happens if the dark en- 
ergy expansion phase preceding the ekpyrotic contraction 
phase lasts at least 600 billion years. Then, a sufficiently 
large patch of space makes it all the way down the classi- 
cal trajectory and through the big bang such that, four- 
teen billion years later, it comprises the overwhelming 
majority of space. This surviving volume, which grows 
in absolute size from cycle to cycle, consists of a smooth, 
flat, expanding space with nearly scale- invariant curva- 
ture perturbations, in accordance with what is observed 
today. As with the mythical phoenix, a new habitable 
universe grows from the ashes of the old. 

The phoenix picture is in some sense minimalist. No 
features have to be added to the cyclic model that were 
not already there, e.g., in the original version inspired by 
M theory and the Hof ava-Witten model 0] . One only has 
had to appreciate the effects of the elements that were in 
place. The heterotic M-theory embedding of the cyclic 
model flj] already admits a four-dimensional effective 
description in terms of gravity coupled to two canoni- 
cally normalized scalar fields. The scalars parameter- 
ize geometrical quantities in the higher-dimensional the- 
ory, namely the distance between the boundary branes 
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and the volume of the internal Calabi-Yau manifold. 
Also, dark energy was always incorporated into the cyclic 
model to stabilize the oscillatory behavior. Here dark en- 
ergy is given the added role of ensuring the survival of 
the phoenix universe. The result is a universe with a 
global structure dramatically different from eternal infla- 
tion and a history somewhat different from earlier oscil- 
latory and cyclic universes. 

For concreteness, we will develop the phoenix picture 
based on the Hofava-Witten model Q so that there is 
simultaneously a 4d and higher dimensional interpreta- 
tion. In this picture, the relevant degrees of freedom aside 
from the two scalars is the radiation (with energy den- 
sity p±) produced on the positive and negative tension 
branes respectively, at the big crunch/big bang transi- 
tion; it appears in the effective theory with a coupling 
(3± that depends on the scalar fields [J] . Thus the equa- 
tions of motion are 

3H 2 = l -^\ + 4,l) + V{&) + P 4 + (fa) P+ + pt(fa) P -{l) 
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3- = + cf>i) + V(fa) - /3+(&W - PL{<l>i)p-, (3) 
a 

where we have assumed a Friedmann-Robertson- Walker 
background with line element ds 2 = — di 2 +a 2 (i)dx 2 ; a(t) 
denoting the scale factor and H = a/ a, with' = Jj. These 
equations form a closed system subject to the equation 
of continuity 



dt 



(P 4 ±P±) + 4Hf3 4 ±P± =Q, 



which implies that 
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Hence we can treat all radiation as one, without the need 
to know the precise functional form of j3± . 

The effective theory contains a potential that is the 
sum of two parts, 



V = V ek (fa,fa) + V rep (fa). 



(6) 



The repulsive potential V rep reflects the fact that there is 
a boundary to scalar field space at fa = 0. This boundary 
corresponds to the negative-tension brane being repelled 
by a naked singularity in the higher-dimensional bulk 
spacetime, and the precise form of the potential depends 
on the nature of the matter living on the negative-tension 
brane (e.g. scalar matter or a perfect fluid) and on its 
coupling to the Calabi-Yau volume modulus; see [13, HU 
for details. A simple form for the repulsive potential is 
V rep oc fa 2 , but we have checked that our results are 
also valid for other functional forms of V rep arising in 
this way. The ekpyrotic potential is the two-field version 
of the cyclic potential described in 
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FIG. 1: This figure shows the evolution of adjacent classical 
trajectories from the ekpyrotic phase through a kinetic phase 
to the big crunch. The trajectories get reflected by a boundary 
in scalar field space at fa = 0. All times are in units of the 
Planck time. 



The essential features are that the potential is steep and 
negative, but with a positive plateau at large field values. 
Here Vo,Vi,Ci ^> 1 are constants, although in general a 
small dependence of the Cj on the fields fa is both natural 
to expect and necessary to obtain a spectrum with a small 
red tilt Q . For the purposes of this paper we can restrict 
ourselves to considering constant Ci . During the ekpyrotic 
phase, the background trajectory is given by the scaling 
solution 



a=(-ty, fa = -ln(-^cfV~/2t), 



(8) 



where t = corresponds to the brane collision and 
p = 2/c 2 + 2/c 2 <C 1 is the equivalent of the slow-roll 
parameter in inflation. The trajectory evolves al ong the 
ridge in the potential. This is most easily seen [lj, [l6| 
by changing variables to a (pointing along the direction 
of the background trajectory) and s (pointing transverse 
to the trajectory): 
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(9) 
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Fluctuations in a correspond to adiabatic perturbations, 
while fluctuations in s are entropy fluctuations. During 
the ekpyrotic phase, the classical background trajectory 
must follow the ridge at s = 0, and around this trajectory 
the potential can be expanded as 

V ek = V -V e-^"(l + -s 2 + ■■■), (10) 

P 

where (by a shift in fa ) we have chosen V\ such that the 
prefactor of the second term is Vq also. This expression 
for the potential clearly shows that the s direction is un- 
stable during ekpyrosis. Hence, unlike for the single-field 
cyclic model, the background solution is not an attractor, 
and it is not clear at first how cycling can remain viable. 
In order to appreciate the severity of the instability, it is 



3 



useful to solve the equation of motion for s in the vicinity 
of the ridge fl3T ] . To leading order in p, it is given by 

s-^s = 0, (11) 

and subject to the boundary conditions s(t e k— b eg) = 
3s /(2t ek -beg), s(tek-beg) = 0, for some constant s , h 
is solved by 

S = s ^j + 2 J 2 ); (12) 

ek — beg 

in other words s grows like the Hubble parameter H. 
Now, the ekpyrotic phase lasts from t e k-beg to t e k-end = 
e~ Nck t e k-beg and, hence, the spread in s grows by a factor 
e N ' k . N e k denotes the number of e-folds of ekpyrosis, and 
in order for the cyclic potential to interpolate between 
the dark energy and the grand unified scales (reaching 
the grand unified scale is necessary for the perturbations 
to have the right amplitude [9]), we need N e k 120. 
Hence the instability causes s to grow by a huge factor of 
about e 120 . Of course this means that at the start of the 
ekpyrotic phase, the field s must be localized close to the 
ridge with great precision. As is already clear from the 
form of the potential (fT!!)) and as previously derived in 
[HI], being in the ekpyrotic phase until t e k- en d requires 
\s(t e k-end) \ ^ P, which translates into 

\s(t ek - beg )\<pe- N '* (13) 

in Planck units. 

After ekpyrosis, the energy density is dominated by 
the scalar field kinetic energies. This kinetic phase lasts 
for about 10 3 Planck times, during which time the trajec- 
tory bends due to the repulsive potential V r e P , see Fig. [T] 
We stop the evolution about one Planck time before the 
crunch, and resume the description at about one Planck 
time after the big bang, see Fig. [2l What we assume 
about the crunch/bang transition is that the scalar field 
velocities are reversed in direction, and that they are in- 
creased in magnitude by a small parameter x, defined as 

a 

4>i,bang = — -\A + X^i,crunch- (14) 

At the brane collision, radiation and matter are pro- 
duced. In fact, if there is slightly more radiation pro- 
duced on the negative-tension brane than on the positive- 
tension one, then the consequence is precisely a slight 
increase in the scalar field kinetic energies, as assumed 
above 0|. (Some forms of scalar couplings to the matter 
fields, 3+, also lead to increased scalar field kinetic en- 
ergies [4j.) This extra kinetic energy helps to overcome 
the Hubble damping due to the radiation created at the 
bang. The constraint that we impose on the radiation 
is that it should not come to dominate until the a field 
has rolled back up the ekpyrotic potential. As shown 
in 0, G3 this translates into a constraint on the reheat 




FIG. 2: After the big bang, adjacent trajectories proceed 
through a second kinetic phase and roll back up the ekpy- 
rotic potential, where the radiation dominated phase starts 
at t ~ 10 15 Planck times. During the ensuing matter and 
dark energy phases, the scalar fields remain almost immobile. 



temperature T r 

T r < (-V(t ek ^ end ))^\ V ° j VW, (15) 

V \Lek — end ) 

which is easy to satisfy. Finally, we note that the radi- 
ation does not contribute to the repulsive potential V rep 
[l5l ]. However, since there can be extra matter produced 
at the bang, which could increase the strength of V rep , 
we have allowed for a fractional increase in V rep of (at 
most) x- 

The extra kinetic energy of the scalars means that the 
a field quickly rolls back up the unstable ekpyrotic poten- 
tial to the (almost-)plateau, and the trajectory deviates 
little from a straight line (note that the evolution is still 
kinetic dominated at this stage). Numerical simulations 
show that during the whole evolution after the ekpyrotic 
phase, the range of values of s (which grows by a fac- 
tor e 120 during the ekpyrotic phase) increases only by 
a factor of about 5. Also, the final range encompasses 
the initial one as long as y/x < 10 -2 , which is a reason- 
able assumption. Numerical simulations also show that 
the radiation and increase in V rep discussed above have 
a very small effect. Hence it is clear now that there is a 
fixed point to the entire evolution over the course of one 
cycle, and this fixed point is located exponentially close 
to the ridge. One might be worried that at the start 
of the dark energy phase, when the field finally turns 
around on the plateau of the potential, the trajectory 
might deviate substantially. This is avoided as long as 
the gradient of the potential in the a direction is much 
larger than in the s direction, i.e. we must ensure that 
|V^(t/Ks| ~ —y/p/2/\s\ » 1 or, in other words, we need 
to have \s(t e k-beg)\ *C y/p. On account of (fT3|) . this is 
clearly satisfied. 

So now consider a small range of s around the fixed 
point: this range grows by about e Nsk over the course 
of one cycle, but, crucially, during this time the uni- 
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FIG. 3: The global structure of the two- field cyclic universe: 
large smooth and flat regions are interspersed with small re- 
gions that have collapsed and have stopped cycling. The tiny 
regions of lighter shade will turn into the entire habitable re- 
gions during the next cycle. 
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/T r ) parameterizes the total growth 



during the kinetic phases and N rac i = ln(T r /To) is the 
number of e-folds of expansion between the start of radi- 
ation domination and today. Hence if 



N de > N e 



3 7 " 



60, 



(16) 



then the gradient of values of s is sufficiently diluted that 
a small initial flat patch will grow into a larger flat patch 
after each cycle. In other words, an initial patch of space, 
with a spread in field values As around the fixed point 
value, will be amplified over the course of one cycle to a 
region containing a larger patch of space with the same 
spread As. Hence, starting from an arbitrarily small ini- 
tial flat patch with the right values of s one can cycle to 
a large, flat and habitable universe. 

In principle, cycling can be eternal to the past and 
future; in this case, the patch of surviving space in a 
phoenix universe would be infinite all along and would 
increase in volume by an exponential factor each cycle. 
An alternative possibility is that the universe had a def- 
inite beginning: for example if it was born in an initial 
quantum creation event, such as in the Hartle-Hawking 
no boundary proposal [18| , or, in the colliding branes pic- 
ture, through the sudden quantum creation from noth- 
ing of a positive- and negative-tension orbifold plane pair 
with random, but smooth, initial conditions. All that is 
required is that there be a finite probability for a patch 
to have the right values of s to grow from cycle to cycle. 
Even if the branes and/or patch begin small, one ends up 
with an arbitrarily large, smooth, flat, habitable universe 
full of matter and radiation as cycles proceed. 

The global structure of the phoenix universe is quite 
different from that of eternal inflation. First of all, 



there is no amplification of rare quantum fluctuations, 
as there is for inflation. The phoenix universe, like gen- 
eral cyclic models, has the property that the smoothing 
phase (dark energy followed by ekpyrotic contraction) 
has a smaller expansion rate, i.e. a smaller Hubble pa- 
rameter H, than the matter and radiation dominated 
phases. Hence, if rare quantum fluctuations keep some 
region in the smoothing phase while typical regions pro- 
ceed to the matter and radiation dominated phases, the 
rare regions lose in volume to the typical ones. The situ- 
ation is reversed with inflation: rare regions kept behind 
in the smoothing (de Sitter) phase are exponentially en- 
hanced and soon dominate the volume. This principle 
underlies eternal inflation. The result is that most of the 
volume in an inflationary universe is empty, uninhabit- 
able and exponentially expanding, punctuated by rare 
bubbles where inflation has ended, some (but not all) of 
which arc habitable. By contrast, most of the volume in 
a cyclic universe is habitable. 

The phoenix universe adds a new wrinkle because of 
the instability associated with the entropic mechanism. 
Namely, only a tiny fraction of the habitable volume in 
any given cycle survives to the next cycle. Most of the 
volume fails to make it through the ekpyrotic smooth- 
ing phase and collapses through a mixmaster crunch Q 
and/or black hole formation (see also |T§|). Nevertheless, 
because the small surviving fraction grows in volume by 
an exponential factor after the bounce (and the failed re- 
gions collapse), most of space during the cycle remains 
habitable. In this sense, habitability is an attractor in the 
phoenix universe. Of course, to some extent the present 
discussion pertains more to metaphysics than to physics 
- the important question is whether the model discussed 
here leads to any distinct observational consequences. 

One observational signature of great promise has been 
discussed recently: due to the very steep ekpyrotic poten- 
tial, ekpyrotic and cyclic models distinguish themselves 
by a substantial and soon measurable amount of non- 
Gaussianity [2(J, Hi], H3| ■ A second potentially observable 
effect is suggested by the above analysis in the case that 
the cyclic universe has a definite beginning: since the 
dark energy phase has only just begun, there should still 
be a significant spread of values of s across the currently 
observable universe. Hence one would expect the dark 
energy not to be equal everywhere in space, but to vary 
by a fractional amount of 

AV _ V ss (As) 2 ^ (As) 2 
V 



V 



p 



(17) 



Since the spread in values of s satisfies As ;$ p and since 
1CP 4 < p < 10~ 2 for models that give observationally 
acceptable levels of spectral tilt and non-Gaussianity, we 
would expect a spatial variation in dark energy of 

^ £ 10- 4 . (18) 

One would expect to be closer to the upper bound if we 
are still in one of the early cycles, and less and less vari- 
ation as we get to older cycles. The largest variation 
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occurs for a flat region that was smaller than the hori- 
zon size in previous cycles and which exceeds the present 
horizon size for the first time during the current cycle. 
However, there is no particular reason to believe that we 
must be in one of the very early cycles. Nevertheless, 
conceptually the spatial variation of dark energy can be 
regarded as a cosmic clock in scenarios where the cyclic 
universe has a beginning. 

Finally, we wish to emphasize the importance of dark 
energy to the phoenix universe. In previous studies of 
the cyclic universe based on a single field where there 
was no instability during the ekpyrotic phase, the role 
of dark energy was confined to ensuring that the cyclic 
solution was a stable attractor solution to the equations 
of motion in the case that the radion overshoots after 
a bounce (or the branes separate too much after a colli- 
sion) . Dark energy is not needed to smooth or flatten the 
universe; this is easily accomplished during the ekpyrotic 
contraction phase alone. To meet the attractor require- 
ment, only two or three e-folds of dark energy domination 
are needed. For the phoenix universe, though, dark en- 



ergy plays a new role: it expands the region of space that 
meets the conditions on s required to complete the un- 
stable ekpyrotic trajectory and rebound as a smooth, flat 
patch with nearly scale-invariant fluctuations in accord 
with what is observed. With too few e-folds of dark en- 
ergy, this patch would shrink from one cycle to the next 
and would not survive. With at least 60 e-folds of accel- 
erated expansion at the current rate (or 600 billion years 
of dark energy phase) , the patch grows from cycle to cy- 
cle, and the phoenix universe is forever reborn. Dark 
energy is, thereby, given a new role and must satisfy a 
new constraint. The new constraint is rather mild quan- 
titatively, but, qualitatively, by controlling the survival 
or termination of the phoenix universe, it may act as a 
selection criterion that may help to explain why the dark 
energy density is so small and yet non-zero today [23j]. 

We would like to thank Jim Peebles and Neil Turok for 
stimulating discussions and comments on the manuscript. 
This work was supported in part by US Department of 
Energy grant DE-FG02-91ER40671. 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 

[7] 
[8] 
[9] 

[10] 
[11] 



G. Lemaitre, Annales Soc. Sci. Brux. Ser. I A 53, 51 
(1933). 

E. Komatsu 



at 



al. 



[WMAP Collaboration], 
larXiv:0803.0547l [astro-ph], 

P. J. Steinhardt and N. Turok, Science 296, 1436 (2002). 
P. J. St einhardt and N. Turok , Phys. Rev. D 65, 126003 
(2002) [arXiv:hep-th/0111098l . 

J. L. Lehners, Phys. Rept. 465, 223 (2008) 



J. Steinhardt 
063514 (2004) 



L. Lehners, Phys. Rept 
|arXiv:0806."T245l [astro-ph]] . 
J. K. Erickson, D. H. Wesley, P 
and N. Turok, Phys . Rev. D 69, 
[arXiv:hep-th/0312009| . 
D. Garfinkle, W. C. Lim, F. Pretorius and P. J. Stein- 
hardt, larXiv:0808.0542l [hep-th] . 

P. Horava and E. Witte n, Nucl. Phys. B 460, 506 (1996) 
|arXiv:hep th/9510209| . 

J. L. Lehners, P. McFadden, N. Turok and 
P. J. Steinhardt, Phy s. Rev. D 76, 103501 (2007) 



arXiv:hep-th/0702153 



P. Creminelli, A. Nicolis and M. Zaldarriaga, Phys. Rev. 
D 71, 06350 5 (2005) |arXiv:hep-th/0411270] . 
R. Durrer, hep-th/0112026 R. Durrer and F. Vernizzi, 
Phys. Rev. D 66, 08 3503 (2002); C Cartier, R. Durrer 
and E.J. Copeland, hep-th/0301198 R. Brandenberger 
and F. Finelli, JHEP 0111, 056 (2001); F. Finelli and 
R. Brandenberger, Phys. Rev. D 65, 103522 (2002); D. 
Lyth, Phys. Lett. B 524, 1 (2002); Phys. Lett. B 526, 
17 (2002); J.C. Hwang, Phys. Rev. D 65, 063514 (2002); 
P. Peter and N. Pinto-Neto, Phys. Rev. D 66, 063509 



[12 

[is; 

[14 

[is; 
[ie; 
[17; 
[is; 

[19 
[20 
[21 
[22 

[23 



(2002); A. J. Tolley, N. Turok and P. J. Steinhardt, Phys. 
Rev. D 69, 106005 (2004); P. L. McFadden, N. Turok and 
P. J. Steinhardt, Phys. Rev. D 76, 104038 (2007). 
A. J. Tolley and D. H Wesley, JCAP 0705, 006 (2007) 
[arXiv:hep-th/0703101| . 

E. I. Buchbinder, J. Khoury and B. A. Ovrut, JHEP 

0711, 076 (2007) [arXiv:0706.3"903l [hep-th]]. 

J. L. Lehners, P. McFadden and N. Turok , Phys. Rev. D 

75, 103510 (2007) larXiv:h ep-th/06 11259l ; Phys. Rev. D 

76, 023501 (2007) arXiv:hep-th/0612026 . 

J. L. L ehners and N. T urok, Phys. Rev. D 77, 023516 
(2008) |arXiv:0708.0743l [hep-th]]. 

K. Koyama and D. W ands, JCAP 0704, 008 (2007) 
[arXiv:hep-th/0703040] . 

J. Khoury, P. J. Steinhardt and N. Turok, Phys. Rev. 

Lett. 92, 031302 (2004) [arXiv:hep-th/0307132] . 

J. B. Hartle and S. W. Hawking, Phys. Rev. D 28, 2960 

(1983). 

J. K. Erickson, S. Gratton, P. J. Steinhardt and N. Turok, 
Phys. Rev. D 75, 123507 (2007) [arXiv:hep-th/0607164j . 
K. Koyama, S. Mizuno, F. Vernizzi and D. Wands, JCAP 
0711, 024 (2007) [arXiv:0708.432T1 [hep-th]]. 
E. I. Buchbinder, J. Khoury and B. A. Ovrut, Phys. Rev. 
Lett. 100, 171302 (2008) [arXiv:0710.5T72l [hep-th]] . 
J. L. Lehners and P. J. Steinhardt, Phys. Rev. D 77, 
063533 (2008) [arXiv:0712.3"779l [hep -th]]; Phys. Rev. D 
78, 023506 (2008) _arXiv:0804.1293 [hep-th]]. 
J.-L. Lehners, P.J. Steinhardt and N. Turok, in prepara- 
tion. 



